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ABSTRACT: Recently the ATLAS collaboration reported a 30 excess in the leptonic- 
Z4+j ets + Ei channel. This may be interpreted in the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM) by gluino pair production with the decay chain g > qqX9 > 
qüZX9, where X9 and X$ denote the lightest and the next-to-lightest neutralinos with 
singlino and bino as their dominant components respectively. After exploring the relevant 
parameter space of the NMSSM by considering the constraints from the ATLAS searches 
for jets + E755 signals, we conclude that the NMSSM is able to explain the excess at 10 
level with the number of the signal events reaching its measured central value in optimal 
cases, and the best explanation comes from a compressed spectrum such as mg ~ 650GeV, 

ei > 565GeV and m; > 465GeV. We also check the consistency of the ATLAS results 
SÉ the null result of jM CMS on-Z search. We find that under the CMS limits at 95% 
C.L., the event number of the ATLAS on-Z signal can still reach 11 in our scenario, which 
is about 1.20 away from the measured central value. 
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1 Introduction 


Since the discovery of a Higgs-like particle in 2012 [1], the search for new physics beyond 
the Standard Model (SM) has been the most important task of the LHC. So far this 
search has covered a wide range of possible signatures of new physics, especially the signals 
of supersymmetry (SUSY), which are characterized by various combinations of missing 
transverse energy En plus jets and/or leptons (electrons or muons). In this direction, 
one impressive result is that the ATLAS collaboration recently reported a 3c excess in the 
channel of two leptons with an invariant mass located around mz, large Eu and at least 
two jets [2]. For 20.3fb^! of integrated luminosity at a center-of-mass energy of 8TeV, the 


number of the observed leptonic-Z events is 29 (summing over electron and muon pairs) 


in comparison with 10.6 + 3.2 events expected in the SM [2]. 

During past weeks, several works appeared to interpret the excess in SUSY [3-7] and all 
of them employed gluino pair production process. As pointed out in [3], in order to enhance 
the signal rate to meet the ATLAS data, the gluino must be relatively light, mg < 1.2TeV, 
and meanwhile have a rich Z-boson yield in its decay. The latter requirement can be 
satisfied in the General Gauge Mediation model (GGM) with a very light gravitino G as 
the lightest supersymmetric particle (LSP) [8]. In this framework, due to the very weak 
couplings of G to matter fields, gluino prefers to decay firstly into a neutralino, and the 
subsequent decay of the neutralino into G is able to produce one or more Z-bosons EI 
The simplest decay chain can be written as à > qqx? > qqZG, where X9 denotes the 
lightest neutralino and the first step decay proceeds via exchanging an off-shell squark. It 
should be noted that, since mg >> mz, one has following mass relations: mg > mgo > mz 
or m3 Z mgo > mz. In either case, the gluino cascade decay will produce at least two 
hard jets, which are from the first step decay g — qqX? or from the hadronic decay of an 
energetic Z boson. As a result, this explanation could be tightly constrained by the search 
for SUSY in events with jets plus E755. Moreover, as illustrated in [4], this explanation can 
not reproduce the shape of the Be distribution measured by the ATLAS collaboration. 

Given the deficiencies of the GGM in explaining the excess, we in this work con- 
sider an alternative simple scenario realized in the Next-to Minimal Supersymmetric Stan- 
dard Model (NMSSM) [10]. Our scenario involves a singlino-dominated LSP X? and a 
bino-dominated next-to-lightest supersymmetric particle (NLSP) X9, and the decay chain 
à — qdX9 — qqZX$. Similar to the GGM explanation, the LSP in our scenario has very 
weak couplings to matter fields, but different from a very light G in GGM (which is re- 
quired by the prompt decay of the NLSP [9]), now the LSP may be massive. This feature 
enables us to choose a slightly compressed SUSY spectrum mg ~ mgg ~ mio + mz to 
evade the constraints from the SUSY searches. Consequently, gluino as light as 700GeV is 
still experimentally allowed, which is helpful to explain the excess (see discussion below). 
Another advantage of our scenario is that it can reproduce the E755 and Hr distributions 
measured by the ATLAS collaboration. This has been recently emphasized in [4]. 

About our scenario, two points should be noted. One is that, since Z boson de- 
cays dominantly hadronically, in most cases the gluino cascade will produce multiple jets. 
Considering that the ATLAS collaboration has performed several comprehensive searches 


for the jets + Eis" signals, and each of the searches focused on multiple signal regions 
which are sensitive to different configurations of SUSY spectrum [17-19], one should con- 
sider these searches by reproducing them through detailed simulations and then using the 
corresponding experimental data as input to limit SUSY. Studying the capability of our 
scenario under these constraints to interpret the excess is the main aim of this work. The 
other is that the CMS collaboration recently also reported its analysis on the leptonic- 
Z + jets + E755 signal [20], which should have the same physical origin as the ATLAS 
excess in our scenario. Since the CMS observed no excess, one should check the consistency 
of the two searches by noting that they are based on different cuts and different detectors. 
'This is another aim of this work. 

This work is organized as follows. In Section II, we introduce our scenario and illustrate 
its key features. Then in Section III, we describe the strategies of the ATLAS and CMS 
searches for the leptonic-Z + jets + Ens signal as well as other searches for jets 4- ae 
signals. In Section IV, we first show the dependence of various signals on SUSY parameters 
in our scenario. Then we scan the relevant SUSY parameters by considering the constraints 
from the SUSY searches, and in the surviving parameter space, we investigate the capability 
of our scenario in explaining the ATLAS Z-peaked excess and check the consistency of the 
measurements from the two collaborations on the leptonic-Z + jets + ges signal. Finally, 
we draw our conclusion in Section V. The validations of our simulations on the ATLAS 


miss 


searches for the leptonic-Z + jets + E^? signal are presented in the appendix. 


2 A potential explanation of the excess in the NMSSM 


'The NMSSM is the simplest singlet extension of the Minimal Supersymmetric Standard 
Model (MSSM) and its superpotential is given by [10] 


AIT ^ ^ l,a Ka 
Wnmssm = Wussm + ASH,- Hat EFS + Sdt + 39^ (2.1) 


where Wmssm is the superpotential of the MSSM, H,, and Hg are MSSM Higgs super- 
fields, is a gauge singlet superfield, A and « are dimensionless Yukawa couplings, W is 
a supersymmetric mass and ép with dimension of squared mass parameterizes a tadpole 
term. 

Due to the addition of the singlet field, the NMSSM predicts three CP-even Higgs 
bosons, two CP-odd Higgs bosons and five neutralinos. The mass matrix of the neutralinos 
in the bases (—iB9, —iW?, H9, H9, S?) is given by[10] 
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View vru 
0 Mə — VOR SC 0 
M = mel SC 0 —Heff —Ava ; (2.2) 
SC SC Hat 0 — Au 
0 0 —Auq —Avy 2KS + p 


where cy, = cos 0w, Mı and Mo are soft gaugino masses, Uu = vsin and vg = v cos B are 
vacuum expectation values (vev) of the fields H. and Hq respectively, and Herr = u + As 


with s denoting the vev of the singlet scalar field S. This matrix indicates that the mixings 
between the bases depend on the parameters tan f and A, and for a singlino-dominated 
LSP, its gaugino components should be very small, while its Higgsino components may be 
sizeable for a large A. 

The interactions of the neutralinos take following form [10] 
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where N is the rotation matrix to diagonalize the neutralino mass matrix with its element 
denoted by Mi; (i,j =1,--- ,5), OF = -O;* = - Nia N} +4 Na ND, and II? = NiaNjo + 
NjyN;,. This Lagrangian indicates that the interactions of the ao with a light 
quark and those with a Z-boson are determined by their gaugino components and Higgsino 
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components respectively, while their interactions with a Higgs boson are decided by the 
both and there may exist cancelation between the two contributions. These features are 
helpful in understanding our scenario. 

As mentioned in last section, the key issues in explaining the leptonic-Z excess are 
how to improve the Z yield rate in gluino decay and how to escape the constraints from 
the ATLAS direct searches for SUSY. Our scheme is as follows: 


e We set the LSP to be the singlino-dominated neutralino, the NLSP to be the bino- 
dominated neutralino, and Mo, uer; > mg. This setup guarantees that the rate of 
the gluino decay into the X9 (via an intermediate off-shell squark) is much larger than 
that of the decay into the X9. Meanwhile, since the NLSP is not accompanied by a 
light chargino, the gluino may decay into the X9 at a rate of approximate 100%. 


e We assume m;,m; > mg > mg where q denotes light quarks. This assumption 
ensures that the gluino decays dominantly into qqX5 state at first step. We suppress 
the decay chain g > HXH — H SX due to its very low efficiency in getting the desired 
signal (We obtain this conclusion by detailed simulations). We assume a relatively 


heavy 4 since it is favored by current LHC search for SUSY. 


As will be shown later, we are particularly interested in the situation mg ~ mx. In 
this case, the assumption on the mass spectrum may be relaxed to be me, mz, mg > ms; 
and mg < mgg + 2m. 


e We require that the decay X9 — Z% proceeds at a rate of approximate 100%. 
'This can be realized by setting mz « mo — men < m or by choosing specific 
SUSY parameters such that the coupling hX9X9 is suppressed. With the package 


NMSSM Tools[11], we numerically checked that, by tuning the parameters of the 


NMSSM, the rate of AN > ZX} can be significantly larger than that of AN — hx? 


even when mao — mxo > Mh. 


e We employ a moderately compressed spectrum mg ^ mgo ~ men + mz to avoid the 
constraints from the LHC searches for SUSY. 


About our scenario, we have more explanations. One is that the measured SM-like 
Higgs boson mass in our scenario can be satisfied by choosing appropriate values of A, tan 8 
and also third generation squark masses [12-15] and the correct dark matter relic density 
can be achieved by the annihilate channels X9X9 — H;/A; — XY with XY being any 
possible final states [16]. The other is that, in our scenario, we ad hoc forbid the possibility 
of M», eff < mg. The main reason is that in that case, the gluino decay chain may become 


more complicated, and channels such as g > qd'X; > qd WA — qq'W-ZX9 are open. 
Then one has to scrutinize the vast parameter space of the NMSSM to find parameter 
points where the Z yield rate is high. Since this process involves the simulation of the 
gluino pair production at the LHC, it is very time-consuming in calculation by cluster 
because the decay chains are lengthy. So our scenario should be regarded as a simplest 
attempt to explain the excess, which only involves the parameters mz, man and mgg. 


3 Simulations and cuts 


In order to discuss the constraints of the LHC searches for SUSY on our scenario and 
also calculate the event number of the leptonic-Z 4-E7"55 signal, we implement the corre- 
sponding experimental searches in the package CheckMATE-1.2.0 [23], which have tuned 
the package Delphes3.0.10 [24] to simulate the behavior of the ATLAS (CMS) detector. 
The validations of our simulations are presented in Appendix for the benchmark points 
provided by the experimental groups, and we find that our results agree with the corre- 
sponding experimental analyses at 2096 level. Throughout this work, the events in our 
simulation are generated by MG5.aMC[25| which includes Pythia|26| for parton shower- 
ing and hadronization, and the package Prospino [27] is used to calculate the NLO cross 
sections of the gluino pair production. 

In the following, we only briefly describe the cut flows of these experiments. Detailed 
information about them can be found from the corresponding experimental reports. 


3.1 ATLAS search for leptonic-Z + jets + Ene signal 


This search concentrated on the events with a same flavor opposite-sign (SFOS) dilepton 
pair, jets and EZ'55 at 8-TeV LHC [2]. One intriguing result of this search comes from 
the signal region SR-Z where the invariant mass of the lepton pair locates around mz. 


Explicitly speaking, in contrast with 4.2 +1.6 electron pair events and 6.4 + 2.2 muon pair 
events expected for the SM background, 16 events and 13 events are observed respectively 
in this region. This implies a 30 excess of the observed events over the background, and 
may be regarded as a hint of SUSY. So in the following we try to explain this excess with 
our scenario. 'l'he signal region SR-Z is defined by 


Table 1: Definition of signal regions in ATLAS search on > 6 jets + E55 signal [19]. 
Here the jet, the pseudo-rapidity 7 and multiplicity all refer to the R = 0.4 jets. Composite 
jets with a larger radius parameter R = 1.0 are used in the multi-jet + M stream when 


constructing the quantity M (see text). 


Multi-jet + flavour stream Multi-jet + M stream 
Name 8j50 9j50 | 210j50 | 7j80 | 28j80 | 28j50 | 29j50 | 21050 
Jet Pr" [GeV] 50 80 50 
Jet count =8 =9 >10 =7 >8 >8 >9 >10 
b-jet count 0,1,>2 | 0,1,>2 0,1,>2 | 0,1,>2 
MŽ [GeV] > 340, > 420 


e The first two leading leptons are SFOS, their transverse momentums are larger 
than 25GeV and 10GeV respectively, and meanwhile their invariant mass satisfies 
81GeV < my < 101GeV. 


e More than one jets have Pr > 35GeV and |n| < 2.5. 


e See > 225GeV and Hr > 600GeV, where Hr is defined as the scalar sum of the 
transverse momenta for all signal jets and the two leading leptons. 


e The azimuthal opening angle between each of the leading two jets and E755 should 
be larger than 0.4, i.e. Ad(jeti2, ERSS) > 0.4. 


3.2 ATLAS search for > 6 jets + En signal 


To search for final states with large jet multiplicities and Eis", ATLAS used two inde- 
pendent streams, i.e. multi-jet + flavor stream and multi-jets+ M> stream, in its anlaysis 
[19]. 

In the multi-jet + flavour stream, the number of jets with |n| < 2 and pr above a 
threshold (50GeV or 80GeV) is firstly determined for a given event. Then events are 
categorized by the number of jets, and further subdivided by the number of jets that are 
b-tagged. After such treatments, totally thirteen signal regions are defined, which are given 
in Table.1. 

Analysis of the multi-jet + M stream seems more complicated. As the first step, the 
number of (R = 0.4) jets with |n| < 2.8 and pr above 50 GeV is determined for a given 
event, and like the former stream, events are classified by their jet numbers. As the second 
step, the four-momenta of the R = 0.4 jets are used as inputs for a second iteration of 
the anti-k, jet algorithm [21], which now adopts a larger jet radius parameter, R = 1.0. 
'The resulting larger objects are called as composite jets. Then a selection variable M is 
defined by 

M? 2 VY mëi, (3.1) 
j 


where the sum is over the mass of the composite jets, and two signal regions are defined 
by the threshold of M (340GeV or 420GeV). In this stream, there are totally six signal 
regions (see Table.1). 

We mention by the way that both the streams veto events containing any isolated 
electron or muon candidates with Pr > 10GeV at the start of the analysis, and require 
des /JHr > AGeV'/2 at the end of the cut flows, where H7 is the scalar sum of the 
transverse momenta of all jets with Pr > 40GeV. 


3.3 ATLAS search for 2 ~ 6 jets + piss signal 


The search for events with 2 ~ 6 jets + E755 at 8-TeV LHC with a total integrated 
luminosity of 20.3fb~! has been carried out in [17] and also in [18]. Compared with the 
preliminary search described in [17], the updated search in [18] defined more signal regions 
to optimize the search. For example, it has been shown that the signal regions 27W and 
4jW are able to improve greatly the sensitivity of the ATLAS search to models predicting 


the cascade g > qq X7 > qd WAN. The main cuts of this search are: 
e E755 > 160GeV. 


e Pr(ji) > 130GeV for the leading jet, and Pr(j) > 60GeV for the other jets. An 
exception is the signal region 4jW where Pr(j) > 40GeV is required. 


e The azimuthal angles between jets and the E755 satisfy óó(jet4 » (3), ER’ ln > 0.4 
and óo(jetisa, Em555),. i. > 0.2. 


e A variety of signal regions is then defined by the number of jets, the values of 
Ess  / Hr, EW /mep¢(Nj) and megr(incl). Here Hr, megg(N;) and me, (incl.) 
are defined as the scalar sum of the transverse momenta for all Pr > 40GeV jets, the 
leading N; jets and E755, and all jets with Pr > 40GeV and E755 respectively. 


In practice, we also consider the search in [17]. This search has been implemented in 
the CheckMATE by the package authors, and we have verified its correctness. Since this 
search is significantly weaker than that of [18] in constraining SUSY parameter space, we 
do not list its cuts here. 


3.4 CMS search for leptonic-Z +jets + E755 signal 


In [20], the CMS collaboration carried out a dedicated on-Z counting experiment to search 
for leptonic-Z +jets + Era signal. This signal is same as that of the ATLAS search in 
[2], but different from the ATLAS result, the CMS collaboration saw no excess. Since the 
two signals have same physical origin in our scenario, they should be correlated. So in our 
work, we consider the CMS search as a possible constraint on our interpretation of the 
excess. The main strategies of the CMS search are given by 


e There exists at least one SFOS dilepton pair (electron or muon pair) with Pr > 
20GeV and |n| < 2.4 for each lepton, and the invariant mass of the lepton pair 
satisfies 81 « my « 101GeV. 


Table 2: Six bins of the CMS dedicated counting experiment for events with an on-shell 
Z boson [20]. For each bin, the 95% CL upper limits on the number of signal events, S9? 


obs? 
are also presented. 


Nyets > 2 23 

EméssiGeV] | 100-200 | 200-300 2300 | 100-200 | 200-300| 2300 
Naata 1187 65 7 490 35 6 
Nokg 12044106 | 74.5411.3 | 12.844.3 | 478443 | 39.246.6 | 5.342.3 
S% 207 20 7.6 89 16.1 8 


e Six bins are defined by the value of E755 and the number of jets with Pr > 40GeV 
and Jul < 3.0. Here we present some information of these bins in TABLE 2. 


About the CMS search, we note that, unlike above ATLAS searches, it does not provide 
the 9596 C.L. upper limits on the number of signal events for each bin, which was denoted 
as S9». in literature [23]. In order to discuss the constraint of such a search, we calculate 
these limits by the asymptotic C'L, prescription [22], and present the corresponding results 
in TABLE 2. As we will show below, the bins with N; > 2,3 and 200 < Emiss < 300 
usually put the strongest limitation on any explanation of the excess, which is different 


from previous analysis [4]. 


4 Numerical results and discussions 


As we introduced in section II, the branching ratio of g > qqX9 > qQZ X4 in our scenario can 
approach 100% by setting the SUSY parameters other than the masses mg, mio and myo at 
certain values. In this case, the event number of the gluino pair production with a certain 


final state is solely determined by the masses, or equivalently by mg, Am; = mg — ma 


xg — mao. In this section, we investigate the constraints of the SUSY 


searches on the masses. For this end, we define for each experimental search the ratio 


and Amo = m 


R = maz(Nsj/ 595%) where Ng, is the signal event number for i signal region in the 


obs,i 
95% 
sobs. 


Obviously, only in the case of R < 1, the corresponding mass spectrum is allowed by the 


search, is its 95% upper limit, and the maz is over all signal regions of the search. 
search at 95% C.L.. We also investigate how well our scenario can explain the leptonic-Z 
excess. After considering the statistical uncertainty of the measured event number, we infer 
from the ATLAS search that the needed signal event number, Ny, is 18.4 + 6.3 (assuming 
the experimental data obey Poisson Distribution and adding the statistical uncertainty and 


systemic uncertainty in quadrature), which corresponds to 12.1 < Ny € 24.7 and 5.8 < 
Ny < 31 in explaining the excess at lo and 2o levels respectively. In following discussion, 
it is convenient to define the normalized leptonic-Z event number as Ry = Nj/18.4, then 
the lo and 2o ranges are restated as 0.66 < Ry < 1.34 and 0.32 < Ry < 1.68 respectively. 
In Fig.l, we show the dependence of Rs for different searches on Am , Am» and mg 
respectively by fixing mg — 4.5 TeV. From this figure, one can learn at least five facts. 
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Figure 1: The dependence of Rs for different searches on Amı, Am 2 and mg. Here 


X X 
that between X$ and X9. For the ATLAS leptonic-Z analysis, R is the normalized lepton 


Amı = mğ — mgg denotes the mass splitting between g and X9 and Ama = mo — mo is 


pair event number, R = Nj/18.4 with 18.4 representing the central value of the excess, 
$9570 


obs,i 


is its 95% upper limit, and the maz is over all signal 


while for any of the other searches, R = max(Ngj/ 


$ H $ . 95% 
for i signal region in the search, SA 


regions in the search. In the latter case, R < 1 means that the parameter point is allowed 


) where Ng; is the event number 


by the corresponding search at 95% C.L.. 


1. The left and middle panels indicate that, as the SUSY spectrum becomes compressed, 
the value of R for each search decreases monotonously. This reflects that a smaller 
mass splitting is able to suppress the cut efficiency in the search so that fewer SUSY 
signal events are finally retained. Consequently, this case is apt to escape the con- 
straints from the direct search for SUSY at the LHC, but meanwhile, the ATLAS 
leptonic-Z events are suppressed too. 


2. Motivated by above observation, we consider a rather compressed spectrum in the 
right panel, and study the dependence of R on mg. The panel indicates that Rs 
decrease monotonously as the gluino becomes heavy. This is because the suppression 
of the gluino pair production rate with the increase of mg usually dominates over the 
enhancement of the cut efficiency in getting the signal events of the search, and the 
net result is the dropping of the signal events. 


3. All the three panels indicate that, in the case of a very compressed spectrum, the 
CMS search for the leptonic-Z +jets + qnss signal puts the severest constraint on 
SUSY among the searches. Especially, as the gluino becomes light, the signal event 
number in the CMS search grew much faster than those of the other searches. As a 
result, the CMS search may be used to set a lower bound on m; in our scenario. 


Miss 


4. In general, the ATLAS preliminary search for 2 ~ 6 jets + ET" signal and that for 
> 6 jets + Ies signal are comparable in constraining the parameter space of our 
scenario, and the both constraints are weaker than that from the ATLAS updated 


search for 2 ~ 6 jets + ERSS signal. 


5. Since the values of Hs for the experimental searches have different behaviors with 
the variation of the mass spectrum, it is better to scan the masses in studying the 
capability of our scenario to explain the excess. 


In Fig.2, we show the constant contours of the event number for the ATLAS search for 
the leptonic-Z +jets + Jue signal on the Am — Am» plane for different gluino masses. 
Regions between line 12.1 and line 24.7 and those between line 5.8 and line 31 are able 
to explain the excess at lo and 2o levels respectively. Constraints from the other SUSY 
searches are also plotted by different types of lines with the left side of the lines being 
experimentally allowed. Solid line, dashed line and dotted line represent the constraints 
from the ATLAS searches for > 6 jets + Eis, 2 ~ 6 jets + ERS (updated) and 2 ~ 6 
jets + E755 (preliminary) respectively, and they are obtained by fixing the corresponding 
R value at 1. The CMS constraint is plotted in a similar way, but by dash-dotted lines. In 
getting this figure, we scan the parameters Am, and Am» with a grid of 30GeV for Am, 
and 10GeV for Amy. 

From Fig.2, one can learn following facts: 


e With the decrease of the gluino mass, only a more compressed spectrum can survive 
the constraints from the LHC searches for SUSY. Among the constraints, the CMS 
search for the letponic-Z +jets + Jupes becomes the most stringent constraint for 
mg S, 800GeV. Especially, in the case of mg = 700GeV, only a very small corner on 
the Am; — Amz plane survives the constraint. 


e If one does not consider the CMS constraint, our scenario can explain well the ATLAS 
leptonic-Z excess without conflicting with the ATLAS searches for jets + Bes. 'The 
central value of the excess (18.4 events) may be achieved for mg < 700GeV, and one 
of the best benchmark points is m; = 650GeV, mgg = 565GeV and mgo = 465GeV. 
While if one take the CMS constraint serious, we find that the leptonic-Z events in 
the ATLAS experiment is at most 11. This tension arises since the CMS signal and 
the ATLAS excess have same physical origin. 
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Figure 2: Constant contours of the event number for the ATLAS search for the leptonic- 
Z +jets + MMs signal on the Am; — Am» plane for different gluino masses. Regions 
between line 12.1 and line 24.7 and those between line 5.8 and line 31 are able to explain 
the excess at lo and 2c levels respectively. Constraints from the other SUSY searches are 
also plotted by different types of lines with the left side of the lines being experimental 


allowed. For the case of mg < 650GeV, there are no regions surviving the constraint from 
the CMS search for the leptonic-Z signal. 
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5 Conclusion 


Recently the ATLAS collaboration reported a 3c excess in the leptonic-Z + jets + EE 
channel. In this work, we intend to explain the excess by the gluino pair production 
process in the NMSSM. For this purpose, we implemented this experimental analysis (and 
also some other useful experimental analyses) in the package Check MATE, and verified the 
reliability of our simulations. 

The scenario we considered involves a singlino-like LSP, a bino-like NLSP and an 
approximate 100% branching ratio for the decay chain à > qqX9 — qqX4. In this case, 
the SUSY signal rates after cuts are essentially determined by the masses mz, mo and 
mo. We scanned these masses by considering the constraints from the ATLAS searches 
jor jets Eier signal, then we concentrated on the leptonic-Z signal in the experimentally 
allowed parameter space. We concluded that our scenario is able to explain the excess at 
Lo level with the number of the signal events reaching its experimental measured central 
value in optimal cases, and the best explanation comes from a compressed spectrum such 
as mg c 650GeV, m. GE c 565GeV and m; BE c 465GeV. We also tested the consistency of 
the ATLAS result Fh the CMS result in the same channel, which observed no excess on 
the leptonic-Z production. We found that under the CMS limits at 9596 C.L., the number 
of the ATLAS signal events can still reach 11 in our scenario, which is about 1.26 away 
from the measured central value. 

We emphasize that, although we employ a specific scenario in our work to explain 
the ATLAS excess, our conclusions seem to be general. That is, the NMSSM is able to 
explain quite well the excess if only the constraints from the ATLAS searches for SUSY 
are considered. While on the other hand, if one considers the CMS search limits, there is 
a mild tension between the ATLAS excess and the CMS result. 

Note added: At the final stage of this work, there appeared a paper by Ulrich 
Ellwanger [4], which considered a similar scenario to ours. Anyhow, our work differs from 
his work mainly in following three aspects: 


e In his paper, only two benchmark points were discussed, while we perform a scan over 
the relevant parameter space. Our calculation by cluster is very time-consuming. 


e We considered more constraints than in his paper, where only the ATLAS preliminary 
search for 2 ~ 6jets + E7755 was considered. Especially, the benchmark point P2 
in his paper was considered as an optimal point to explain the excess, while in our 
work, we got R = 1.19 for the updated ATLAS search for 2 ~ 6jets + Eier (with 
mg = 4.5TeV). Considering that we have about 20% uncertainty in simulating 
thesearch, we infer that the point is at the boundary of being excluded. 


e When we consider the CMS analysis, we found that the most stringent constraint 
comes from the bins of Njets Z 2,3 and 200GeV < Eee < 300GeV for most cases, 
instead of the bin Njets 2 2 and des > 300GeV adopted in [4]. This is because 
we are considering a compressed spectrum, where the press of the signal is usually 
small. 
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A Validation of Monte Carlo simulations 


In this appendix, we show the validations of our simulations on the ATLAS searches for 
the leptonic-Z +jets + E755 signal [2] and the 2 ~ 6 jets + EZ"55 signal [18]. For each 
search, we choose two benchmark points adopted by the ATLAS report, and repeat the 
experimental analysis. Our simulation results are presented in Table 3 and Table 4 respec- 
tively, and all of them indicate that our simulation agrees with the corresponding ATLAS 
analysis within 20% uncertainty. We also validated our simulation on the CMS search for 
the leptonic-Z +jets + Eier signal [20] with the corresponding results given in Table 5. 
This table indicates that we can reproduce the CMS results within 1096 uncertainty. 
Since the ATLAS preliminary search for 2 ~ 6 jets + E755 signal [17] and the search 
for the > 6 jets + E755 signal [19] do not affect our conclusions, we do not present their 
validation here. In fact, we checked that we can reproduce these searches within 1596 


uncertainties. 
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Table 4: Validation table of our implement of the ATLAS search for the 2 ~ 6 jets + 
E755 sional [18] in CheckMATE. We considered the process pp > 99 > Xi X13Jjj — 
XOX0W-*Wjjjj at 8-TeV LHC in the MSSM with two benchmark points taken from [18]. 
Numbers in the columns EXP and OUR are the cut efficiencies in 27W or 6jt signal region 
obtained by the ATLAS collaboration and us respectively, and those in the column DIFF 
represent their relative difference. 


g g one step, mg = 1200GeV, Mot = 1150GeV, mgo = 60GeV 
SR:2;W EXP | OUR | DIFF 
E555 > 160GeV, Pr(ji2) > 130(60)GeV | 52.70 | 56.35 | -7% 
Ad¢(j1,2,3, Emis) > 0.4 46.30 | 49.01 | -6% 
N(W) unresolved > 2 9.20 | 8.70 5% 
ERs ie, ¢ ¢(Nj) > 0.25 7.00 | 6.69 4% 
Mer p (incl.) > 1800GeV 5.30 | 4.86 8% 

g g one step, mg = 1265GeV, mgt = 945GeV, mo = 625GeV 
SR:6jt EXP | OUR | DIFF 
E755 > 160GeV, Pr(ji,j2) > 130(60)GeV | 53.30 | 54.16 20 
Pr(j3) > 60GeV 53.00 | 53.83 | -2% 
Pr(j4) > 60GeV 50.50 | 51.50 | -2% 
Pr(js) > 60GeV 41.40 | 43.49 | -5% 
Pr(je) > 60GeV 26.70 | 29.78 | -12% 
Aó(j1,2,3, ERSS) > 0.4 22.40 | 25.38 | -13% 
Aó(jisa, ER?) > 0.2 18.20 | 20.54 | -1396 
E555 lm, ¢(Nj) > 0.25 10.90 | 11.64 | -796 
Meg f (incl.) > 1500GeV 4.20 | 4.81 | -15% 
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Table 5: Validation table of our implement of the CMS search for the leptonic-Z +jets + 
E755 signal [20] in CheckMATE. We considered the process pp > 99 > XQXijjjj — 
GGZZjjjj at 8-TeV LHC in the GMSB model with the benchmark point taken from [20]. 
Numbers in the columns EXP and OUR are the event numbers of different bins obtained 
by the CMS collaboration and us respectively, and those in the column DIFF represent 
their relative difference. 


EXP OUR DIFF 
All events 37.7 37.7 0 
> 2 leptons(I*/*),Pr > 20GeV 11.9 11.8 -0.5% 
81 < my < 101GeV 10.7 10.4 -2.7% 
Tets > 2 Niets > 3 
EXP | OUR | DIFF | EXP | OUR | DIFF 
Njets > 2 or 3 10.7 | 10.4 | -2.9% | 10.4] 10.2 | -1.9% 
Em ss > 100GeV 10.3 | 10.0 | -3.396 10| 9.8 | -2.1% 
Emss > 200GeV 9.2 | 88 | -4.1% |} 89) 8.7 | -2.5% 
Em ss > 300GeV 7.6) 7.2] -5.4% | 74] 7.1 | -4.3% 
100GeV < ER's < 200GeV 11| 11 3% || 11] 1.1] 0.8% 
200GeV < Eier < 300GeV 16] 1.6 2% || 15] 16] 6396 
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